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ABSTRACT 

The spectral energy distributions (SED) of dusty galaxies at intermediate redshift may look similar to very high redshift galaxies in 
the optical/near infrared (NIR) domain. This can lead to the contamination of high redshift galaxy searches based on broad band 
optical/NIR photometry by lower redshift dusty galaxies as both kind of galaxies cannot be distinguished. The contamination rate 
could be as high as 50%. This work shows how the far infrared (FIR) domain can help to recognize likely low-z interlopers in an 
optical/NIR search for high-z galaxies. We analyse the FIR SEDs of two galaxies proposed as very high redshift (z > 7) dropout 
candidates based on deep Hawk-I/VLT observations. The FIR SEDs are sampled with PACS/Herschel at 100 and 160 jim, with 
SPIRE/Herschel at 250, 350 and 500 yum and with LABOCA/APEX at 870yum. We find that redshifts > 7 would imply extreme 
FIR SEDs (with dust temperatures > 100 K and FIR luminosities > 10 13 L Q ). At z~2, instead, the SEDs of both sources would be 
compatible with that of typical ULIRGs/SMGs. Considering all the data available for these sources from visible to FIR we re-estimate 
the redshifts and we find z ~ 1.6-2. 5. Due to the strong spectral breaks observed in these galaxies, standard templates from the 
literature fail to reproduce the visible-near IR part of the SEDs even when additional extinction is included. These sources resemble 
strongly dust obscured galaxies selected in Spitzer observations with extreme visible-to-FIR colors, and the galaxy GN10 at z = 4. 
Galaxies with similar SEDs could contaminate other high redshift surveys. 

Key words. Submillimeter: galaxies - Galaxies: high-redshift - Galaxies: distances and redshifts 



1. Introduction 

Observing galaxies up to very high redshifts allows us to 
study directly the formation and evolution of structures in 
the expanding Universe. Finding galaxies at ever higher red- 
shifts has therefore become one of the main areas of extra- 
galactic astronomy. The most common technique is to use 
known broad features in the spectral energy distributions (SEDs) 
of galaxies to identify high redshift sources in deep optical 
and near infrared (NIR) multi-band observations. In particu- 
lar the Lyman break is widely used to select sources by red- 
shift noting their disappearance in bands below a given wave- 



length, the so-called drop-out technique (Stei del et alj H996). 
With this technique and state of the art telescopes and instru- 
ments it is now possible to select sources that are good can- 
didat es for being at the end or within the epoch of reioniza- 
tion dRichard et alJ2006tlZheng et alJ2009HMcLure et al.ll20irj 
IWilkins et al.ll20ld:IOesch et al.ll20ldlBouwens et al.ll2010blah . 

Low redshift galaxies, however, can have very steep SEDs 
resembling a break in the UV/optical/NIR. This can lead to con- 
tamination of the dropout selection of very high-z galaxies, and 
hence to erroneous estimates of the star formation rate den- 
sity, stellar masses, and others, although these effects are cur- 
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rently difficult to quantify. Such object s have been found and 
discussed by several authors (see e.g. iDickinson et al.l [2000: 
Mob asher etaH 120051: ISchaerer et all 120071: iDunlop et al.1 120071 
Chary et alJ2007tlCapak et alJ201ll) ! Confirming the photomet- 



ric redshifts of high-z galaxies by identifying spectral lines is 
challenging because the sources are generally too faint for spec- 
troscopic follow up observation s or because they m ay intrinsi- 
cally lack Lya emission (but see lVanzella et alJ uOlO). 

The recent developments of space far infrared instrumenta- 
tion offer new perspectives in this domain. In particular, with 
the advent of the Herschel Space Observatory it is now pos- 
sible to sample the Far Infrared (FIR) SEDs of the galaxies, 
where the thermal dust emission dominates. The shape of the 
FIR SED universally looks like a broad bump which can be 
used to further constrain the optical/NIR photometric redshifts. 
Although the wavelength of the FIR SED peak also depends on 
the dust temperature, the limited range of average temperatures 
observed so far in galaxies ( between 20 and 60 K averaged over 
the entire galaxies, see, e g..lKovacs et al. 2006; Magnell i et all 
l20iri fMagdis et al. 2010: Wardlow et alJ2011h can be used as a 
prior and makes it possible to discriminate between intermediate 
(z < 3) and very high redshifts (z > 6). 

FIR observations of high-z candidates are also essential to 
characterize their star forming and dust properties and thus in- 
terpret correctly their contribution to the cosmic history of star 
formation an d reionization . 

Recently Lapo rte et alJ d201 ll) identified ten z > 7 candidates 
in the field of the cluster Abell 2667 using photometric dropout 
criteria based on deep observations with HAWK-I on the ESO 
Very Large Telescope (VLT). Comparing to other studies and 
in particular to the WIRCAM Ultra Deep Field Survey (WUDS; 
Pello et al in prep), which is based on deeper optical observa- 
tions bluewards to the I-band, they estimated that 50-75% of 
these candidates could in fact be lower redshift interlopers. Here, 
we study two galaxies of this sample that are clearly detected by 
Herschel, namely the sources named 'zl' and "Y5\ Our goal 
is to determine whether they could be interlopers and to under- 
stand their nature. The redshift probability d istributions of these 
two sources derived bv lLaporte et al.1 (1201 ll) from SED fitting to 
deep Optical/NIR photometry show a prominent peak at z =7.6 
and 8.6 respectively. However, a secondary peak at lower red- 
shift around z ~ 2, i ndicates that they also could be interlopers. 
Lapo rte et al.l d201 ll) also noted that the 24 fim detection of zl 
with MlPS/Spitzer (Y5 is out of the Spitzer map) seems diffi- 
cult to reconcile with the high-z solution. We use new Herschel 
and LABOCA observations of Abell 2667, to reconstruct the FIR 
part of their SEDs. We can thus further constrain their redshifts 
and study their physical properties. 

The layout of the article is as follows: Section|2] gives a pre- 
sentation of the observations and data analysis. In Section[3] the 
FIR part of the SEDs is analyzed. In Section|4] the complete 
SEDs are used to estimate the redshifts and discuss the phys- 
ical properties of the two galaxies. In Section|5] we compare 
the two galaxies to other similar galaxies found in the literature. 
Section|6]gives the conclusions. We assume a A-cosmology with 
H = 70kms-' Mpc" 1 , Q. M = 0.3 and Q A = 0.7. 



2. Observations and Data Analysis 

2.1. Observations and reduction 

Herschel observations were obtained in the fram ework of the 
Hersch el Lensing Survery (HLS) described by lEgami et al.l 
d2010l) . They include PACS data at 100 and 160/rni, and SPIRE 



observations at 250, 350 and 500 fim. The d ata reduction was 
don e with the HIPE so ftware as described bv lRex etal ] (120101) 
and lRawle etal] d2010h . 

The large APEX Bolometer Camera (LABOCA 

ISiringo etal] £009) is a bolometer array operating at 870 jum 
and m ounted on the APEX telescope in the desert of Atacama, 
Chile dGiisten et all 120061) . The LABOCA observations were 
conducted during the summer 2010. The cluster was mapped 
in spiral mode during 30 hours covering a circular field of ~6' 
in radius. The data were reduced with the BoA[3 software. The 
noise is not uniform over the map and the RMS is in the range 

I. l-3.0mJy, the highest values being reached at the edges of the 
map. 

We also obtained a VLA 1 .4 GHz continuum map of Abell 
2667 (PI: R. Ivison) with an RMS of 46 /Jy. 

2.2. Analysis 

The astrometry of all the maps was corrected to align them 
wi th the VLT Ks ima ge. All the optical drop-out sources 
of [Laporte et al] d201 lb were inspected in the Herschel and 
LABOCA images. Two of them, zl and Y5, are detected in 
several FIR bands. As IRAC/Spitzer and MlPS/Spitzer data are 
available for zl, the source can be followed from one band to the 
next one by increasing wavelength despite the decreasing reso- 
lution. Its identification is therefore robust. 

For Y5 there is a larger gap in the SED due to the lack of 
data between 8 fim and lOOjt/m, and due to the fact that it lies at 
the noisy edges of the 100 and 160 fim maps, where it is not de- 
tected. However, Y5 is the only source detected at 4.5 fim within 
a radius of 3" (i.e., ~ 1 /3 rd of the 250 fim beam radius) around 
the 250 fim peak, its identification with the SPIRE detection is 
therefore very likely. 

The fluxes are measured at the positions of the two galax- 
ies by PSF fitting in apertures with a radius equal to FWHM/3, 
where FWHM is the PSF full width at half maximum, i.e., 5.6", 

II. 3", 18.1", 24.9", 36.6" and 22.5" from 100 to 870jt/m. The 
last (LABOCA) FWHM corresponds to the APEX beam con- 
volved by a gaussian of 12". The sources were deblended from 
the neighbouring sources by subtracting PSFs at the positions of 
the neighbours derived from the 250 fim map. Observations at 
these wavelengths with these resolutions are affected by source 
confusion. As a consequence a measured flux cannot be di- 
rectly interpreted as the true flux of a single underlying source. 
A correct treatment of the effect of source confusion on flux 
measurements (a.k.a. flux 'deboosting') requires a prior knowl- 
edge of the source counts towards low fluxes at the given wave- 
length . We followed the method presented by Crawford et al. 
(12010 ) based on a Bayesian analysis. For the prior sou rce counts 
we ex trapolated toward low fluxes the res ults of iBerta et al] 
J2010h for PACS bands , lOliver et alJ d20~iol) for SPIRE bands 
and lCoppin et al.1 d2006l) for the LABOCA band. 

Blending affects zl photometry at A > 250 fim and Y5 pho- 
tometry at A > 500 fim. And the effect of deboosting is small 
(< 20%), except for the 870 fim measurement of Y5 that corre- 
sponds to a 2.6<r signal and that we chose to consider as a tenta- 
tive detection. The deboosted flux of Y5 at 870yum is 1.8+1 mJy 
for a measured flux of 2.5+0.95 mJy. 

None of the two sources are detected in the VLA map. 
The measured FIR fluxes of th e sources as well as their opti- 
cal to near-IR photometry from lLaporte et al] d201 lb are listed 
in TableQ] Postage stamps of the Herschel and LABOCA bands 



1 http : //www . apex- telescope . org/bolometer/laboca/boa/ 
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Table 1. Multi-wavelength SED of zl (col 3) and Y5 (col 
4). Optical and IR photometry (rows 1-10) is taken from 
lLaporte et al.l (1201 ll) . Upper limits are 3 cr. Rows 12-17 give the 
Herschel and LABOCA source flux. When the measured flux 
(not the estimated flux) is < 2.5 cr the 3<x value is given as an 
upper limit. No entry indicates the lack of data. All fluxes are 
given in milli-Jansky. 
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centered at the source positions as well as FIR SED fits are 
shown in the Fig.Q] Y5 is close to the border of the PACS maps 
where the noise is higher, hence the high upper limits. 



3. Analysis of the FIR SEDs 

The following models or templates have been fitted to the FIR 
measurements (cf. Fig.[TJ): 

- A modified black body SED parameterized as described by 
iBlain et all d2003h . with emissivity fixed to /} = 1.5 and 
the Wien correction parameter a = 2.9. T hese values are 
adapt ed to submillimeter galaxies (SMGs; IChapman et alJ 
20051) and local u l traluminous infrar ed galaxies (ULIRGs; 



Dunne et 

IDEOOO; 

IBlain et aT]|2003l) . The free parameters 



are the total FIR luminosity, Lfir, defined as the luminosity 
emitted in the range 8-1000 yum, and the dust temperature, 
T d . 

- The 105 galaxy templates built by Chary & Elbal (1200 lh . 
The templates are fitted wi thout resca ling. 3 

- ULIRG templates built bv lVega et all (12008). with a scaling 
parameter, Lfir. 

- The starb urst, Seyfert, a n d activ e galactic nuclei (AGN) tem- 
plates of iPolletta et alJ d2007l) . with a scaling parameter, 

- Templates built by [Michalows ki et alJ (l2010allbl) to fit high 
redshift galaxies with detected but poorly sampled submm 
emission, with a scaling parameter, Lfir- 

- The SED fit to the observations of SMM J2135-0102 
dSwinbank et aT]|20 1 Ob [Tvison et al. 2010) with a scaling pa- 
rameter, Lfir- 

The fit is performed by finding the maximum likelihood as- 
suming gaussian probability distributions for the measurements. 
When there is no detection the 3-cr value is used as a hard upper 



limit, i.e., the probability is assumed to be uniform in the [0, 3<x] 
interval an d zero outside . The re dshifts are fixed to the solutions 
derived by lLaporte et al.l (1201 lh from the optical/NIR photome- 
try, i.e., z =1.8 and 7.6 for zl and z =1.7 and 8.6 for Y5. The 
corresponding magnification factors are /u =1.12 and 1.17 for z 1 
and ju =1.04 and 1.15 for Y5. The MlPS/Spitzer 24fim flux of 
zl (Y5 has no 24 fim data available) was taken into account to 
fit various galaxy templates, but ignored to fit the modified black 
body since it is most likely dominated by polycyclic aromatic 
hydrocarbons (PAHs). 

For both sources we find reasonable fits at low redshift for 
the modified black body a nd the various ga laxy templates. A 
ULIRG template from the IVega et aTl d2008l) library, a submil- 
limeter detected galaxy template from the Michal owski et al 
(2010a b) library and a Seyfert template from the dPolletta et al 



2007) library are able to reproduce the 24jt/m emission of zl. 
The modified black body model gives dust temperatures of 34 
and 40 K, for zl and Y5 respectively, which are typical val- 
ues for integrated dust temperatures in LIRGs. The infrared 
luminosities, Lfir, are in the range (1.2 - 1.7) x 10 12 L G and 
(3.1 -4.7) x 10 12 L for zl and Y5 respectively. There is a notice- 
able agreement in Lfir between the modified black body model 
and the various templates. These galaxies would therefore be 
typical ULIRGs/SMGs at z ~ 2. This is consistent with the gen- 
eral picture of galaxy evolution now widely observed, i.e., that 
the contribution of ULIRGs to the cosmic SFR is expected to 
peak at z ~ 2 where it sho uld be comparable to that of the more 
'normal' galaxies (see e.g. jMurphy et aljl201lb . 

For the high redshift solutions (z > 7.5), instead, the mod- 
ified black body requires for both sources very high dust tem- 
peratures, i.e., 105 and 129 K for zl and Y5, respectively. Such 
high temperatures averaged over an entire galaxy are extreme. 0. 
This can be seen from the impossibility to find any good fit in 
the different template libraries, which were built from observed 
galaxies. However, the dust properties of galaxies at such high 
redshifts are unknown and dust temperatures above 100 K cannot 
be ruled out. This would imply that the FIR luminosities are of 
the order 0.5 and 1.5 xlO 14 L for zl and Y5, respectively, i.e. 
both sources would be classified as Hyper Luminous Infrared 
Galaxies (HyLIR G). While the nat ure of HyLIRGs is still a mat- 
ter of debate (e.g. iRuiz et al.l2010l) and their density at very high 
redshift is not well known, they are extreme sources with a lower 
number density than ULIRGs. 

The radio continuum upper limits are too high to constrain 
the SED fitting. We find that for Y5 only at the low redshift solu- 
tion (z ~ 2) the Chary & Elbaz template is close to the 3<x upper 
limit. 

Thus, in summary, by comparing the FIR photometry to 
known galaxy SEDs and by taking into account the expected 
temperature and luminosity range of high redshift galaxies, the 
very high redshift solutions derived from the optical/NIR pho- 
tometry seem to be less likely than the low redshift solution. 
The two sources are most likely typical ULIRGs at z ~ 2. This 
result puts strong constraints on the optical/NIR analysis which 
gave a much higher probability to the very high redshift solu- 
tion when no prior luminosity function was taken into account. 
The FIR data alone, however, cannot be used to derive any ac- 



2 We found that such an additional scaling parameter was not re- 
quired to obtain good fits. 



3 To our knowledge there was only one such extreme ca se reported 
so far: the host of the lense d quasar APM08279 at z = 3.9 (Wei fi et al.1 
l2007l : lRiechers e t al. 2009) . It requires two dust components at 75 and 
220 K and the highest temperatures most likely result from heating by 
the powerful quasar 
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Fig. 1. The top figures show the postage stamps in the five Herschel bands (100, 160, 250, 350, 500 //m from left to right) and the 
LABOCA band (870 /jm, the right-most stamp) centered on zl and the best fit SEDs for the low redshift (left) and high redshift 
(right) solutions. The bottom figures show the same for the source Y5. The red circle in the postage stamps is meant to guide the 
eye, it has a radius of 27" and it is centered on the source position. The blue curves correspond to the best fit modified black body 
SED s with the pa rameters written in blue, the magenta curves correspond to the SM MJ2 135-0 102 m odel dSwinbank et alj 2010; 
llvison etal I l2010h . the ot her curves correspond to th e best fit templates of the Chary & Elbaz (2001) library (red), the IVega et al.1 
(2008j) library (green), the Micha towski et al.l d2010alfbT) library (cyan) and the Polletta et al.l d2007l) library (orange). The names of 
the best fit templates as well as Lfir are written on the figure with the same color codes. 



curate photometric redshift because of the redshift-temperature 
degeneracy. 



4. Analysis of the complete SED from visible to FIR 

We will now examine all the data from the visible to the FIR 
ranges to improve the redshift estimate of our galaxies and to 
examine the nature and physical properties of these sources. 
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4.1. Method 

To model the SED of the two sources we use a modified ver- 
sion of the Hyperz phot ometric redshift code o f lBolzonella et al.l 
(l2000h described in ISchaerer & de Barrosl (I2009I) . Non- 
detections are treated as the usual case 1 of Hyperz, i.e. the flux 
in these filters is set to zero, with an error bar corresponding 
to the flux at lcr level. The basic sp ectral templates are take n 
from the Bruzual & Chariot models (iBruzual & Charlotl l2003). 
computed for a variety of star-formation histories and metallici- 
ties. Although applicable only to a limited part of the spectrum, 
we use these templates here to constrain redshift, extinction, and 
stellar mass in particular. For the Bruzual & Chariot templates 
we consider variable ex tinction with Ay up to 8 magnitudes for 
the|Calzetti et al. (2000) attenuation law. We have also explored 
other extinction laws. 

The code, initially designed to fit rest frame UV to near-IR 
(stellar) emission, can also easily be used to include the thermal 
mid-IR and beyond. To cover the entire spectral range from the 
visible to the millimeter domain, and to compare our sources 
with SEDs of very different galaxy types, we have compiled 
a larg e variety of spect ral templates fr om the GRASIL mod - 
els oflSilva et al.l (fl998). the library of ICharv & Elbazl (l200lh . 
Rieke et al.l (120091). the starburst, Seyfert, and AGN templates o f 
Polletta etal.ld2007l) . the ULIRG te mplates oflVega et al.l (12 008). 
the sub-mm galaxy templates of Michalowski et al.|~(2010a ) 



and the model fit to SMM J2135-0102 (Swinbanket al. 2010; 
llvison et"aTll2010h . Extinction can also be added to these spec- 
tral templates; SED fits with and without additional extinction 
will be discussed below. 

We have carried out both fits of the entire SED (optical, near- 
IR, and IR) and fits up to 8 /mi only (for the Bruzual & Chariot 
templates). For each template set the free parameters are redshift 
and (additional) Ay. Physical parameters such as the infrared lu- 
minosity, Lir, defined as the luminosity emitted in the range 2- 
1000 fim; the IR star-formation rate, SFR; and the stellar mass, 
are subsequently derived from the best-fit templates. In contrast 
to the IR fits discussed in Sect. [3] we have no handle on the dust 
temperature, since this is not a parameter describing the SEDs 
used here. We have also checked that the two independent fitting 
methods used here and in Sect.[3]give consistent results. 

4.2. Photometric redshifts 



As discussed in depth bv lLaporte et alJ d201 ll) the best-fit photo- 
metric redshifts of our sources derived from the optical to near- 
IR photometry (up to 8 fim) and using standard spectral tem- 
plates is consistently found at z > 7 with a lower probability at 
low z. This result remains unchanged with the exploration of a 
wider range of extinction, different attenua tion/extinction laws , 
and templates sets used here compared to Laporte et al] d201 ll) 
On the other hand, analysis of the IR SED and other arguments 
clearly favour low redshifts (z ~ 1.5-2.5) as discussed above. 
Subsequently we therefore limit ourselves to z < 4 and attempt 
to refine the photometric redshift of the two sources. 

4.3. Results forzl 

Overall the global, multi-wavelength SED fits for this source 
are rather satisfactory, as shown in Fig. [2] albeit there are im- 
portant discrepancies in the optical dom ain (cf. below). Se veral 
templates (i.e. Polletta's, the U LIRGs of IVega et all (I2008I) . and 
the best-fit SMG template fromlMichalows ki et al.l(l2010al) ) also 
reproduce the 24 /mi flux, and the observed 100 and 870 fim 



fluxes are within 2-3 cr of the model. Interestingly, the best fits 
for both sets of Polletta's and Michalowski's templates are found 
with templates for active galaxies. 

The best-fit redshift found with these templates are between 
z ~ 2.24 and 2.57. The resulting IR luminosity is Lir ~ (2.6 - 
3.2) x 10 12 L H , the correspond ing SFR * 450-550 M yr _1 using 
the standard lKennicutll (ll 998 ) calibration. Fits to the IR part with 
the SMMJ2 135-0 102 template yield z p h n t ~ 2.0. A somewha t 
lower redshift of z « 1 .7 is found with Bruz ual & Charlotl (120031) 
templates using the SED up to the IRAC bands. The estimated 
extinction is A v ~ 2.6, the stellar mass M+ ~ 6 x 1 10 M for 
the same Salpeter IMF adopted by iKennicuttl ([1998). However, 
these values should be taken with caution as the fits are not of 
good quality. For comparison from the absolute H-band magni- 
tude (M B « -23.0) one obtai ns M+ ~ 3 x 1 10 M using the 
mass-to-light ratio adopted bv lWardlow et al.l (1201 ll) for SMGs. 

At a more detailed level (see right panel), all spectral tem- 
plates have some difficulty to reproduce the steep, observed SED 
between the visible (I, z bands) and the near-IR (Y and J here), 
and they predict a flux excess in the optical domain. Below we 
will show that this also holds when variable extinction is added 
to the empirical templates. The same is also true for all other 
templates we have examined, in cluding the theoretical galaxy 
templates of IBruzual & Charlotl |2003). This sharp drop is of 
course the reason why this source was selected as an optical 
dropout (Y-drop). 



4.4. Results for Y5 

For this source the global fits are less good than for zl. This 
is due to the fact that Y5 shows a higher flux ratio between the 
thermal-IR and the near-IR than zl, whose SED already required 
templates with extr eme IR/near-IR fluxe s. For example, the Arp 
220 template from [Polletta et alJ (120071) . shown in red, under- 
predicts the IR flux by a factor ^ 5. The only template com- 
ing near the observed IR emission is from the SMG library of 
iMichalowskietail d2010al) (SMMJ221725.97+001238). With a 
best-fit redshift of z p h ot ~ 2.15 this translates to Lir = 2.2 x 10 12 
L , corresponding to SFR a 380 M yr _1 . Fits to the IR part 
with the SMM J2135-0102 template yield z ph ot ~ 1.8. A best- 
fit redshift of z * 1.95 is found with IBruzual & Charlotl d2003l) 
templates using the SED up to the IRAC bands. The estimated 
extinction is A v ~ 1.6, the stellar mass M* ~ 3 x 10 10 M . 
However, these values should be taken with caution as the fits in 
the domain close to the optical are not of good quality. Again, 
using the absolute H-band magnitude (Mb ~ -22.3) one obtains 
M^ ~ 2 x 10 10 M with the assumptions already mentioned 
above. 

As for zl, the visible-near-IR drop of the SED (see right 
panel) is not well fit by the spectral templates, predicting that the 
source should be detectable in the visible (I, z, Y bands in par- 
ticular), in contrast to our observations. The templ ate fitting best 
this pa rt of the spectrum is an SO template from Polletta et al. 
(2007), shown in green. However, this template underpredicts 
the IR emission by several orders of magnitudes. 



4.5. Possible explanations for the strong SED break/very red 
spectrum 

As already seen, the common, observed spectral templates fail 
to reproduce the steep, observed SED between the visible (I, 
z bands) and the near-IR (Y and J here), and they predict a 
flux excess in the optical domain. What causes the sharp ob- 
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Fig. 2. Fits to the observed SEP o f source zl (photometry shown by blue symbols, including 3 cr upper limits) using dif ferent spectral 
templa tes: ICharv & Elbazl (12001 1) (black = best-fit template, and magenta = template with maximum IR lumino sity! iPolletta et al. 
(2007) (red — global best-fit template, and green = best-fit to visible-near-IR SED excluding the thermal IR), and Michalo wski et al. 
(2010a) (blue). The best-fit SED with the templates of Vega et al. (2008), very similar to the one using Polletta's templates, is not 
shown here for simplicity. Left: Global visible to sub-mm SED. Right: Zoom on vis ible to near-IR part of the SED including for 
comparison also the best-fit SED at high redshift (z = 7.5) from lLaporte et al] (1201 ll dashed line), which is most likely excluded 
due to our Herschel detections. 

A2667 - Y5 A2667 - Y5 




log(wavelength) [^m] log(wavelength) [Mm] 

Fig. 3. Same as Fig. 0for Y5. The best-fit SED at high redshift (black dashed line) is for z = 8.28. 



served decrease of the flux between the near-IR and the optical 
for these sources? The main difficulty provides from the fact 
that the largest spectral break known in galaxy spectra is the 
Lyman break, whereas the typically observed Balmer (or "4000 
A") break is smaller than that of our two galaxies. We have ex- 
amined various possibilities, but with no convincing answer. 



For example, as shown in Fig. [4] for zl, adding variable ex- 
tinction to the empirical templates allows one to diminish some- 
what their optical flux excess. For the Chary & Elbaz templates, 
the best-fit (to the domain shown here) is then with an addi- 
tional extinction AAy = 1.4 for the Calzetti law. However, the 
templates are too smooth to reproduce the apparent break. The 
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A2667 - zl 




log(wavelength) [jum] 

Fig. 4. Observed visible to near-IR part of the SED of zl 
(blue circles) and best-f it models to this spectral range: us- 
ing iPolletta et al.l (2007) templates and no additional ex tinc- 
tion (green line), b l ack li nes: using ICharv & Elbazl j200Th and 
Bruzual & Chariot (2003) templates with variable attenuation 
described by the Calzetti law, magenta: same as black lines, but 
adoption the SMC extinction la w, red line (continuum) an d red 
crosses (synthetic flux in filters): iBruzual & Chariot! (12003) tem- 
plates with nebular emission and SMC law, black dash ed: high 
redshift (z = 7.5) template from lLaporte et all (1201 ll dashed 
line). In all cases the best-fit redshift is z ~ 1.5-1.7, except for 
Polletta's templates. 



flux at optical wavelengths can further be reduced by assuming a 
steeper attenuation law than the Calzetti law adopted by default. 
With the SMC extinction law bv lPrevot et all dl984h the best-fit 
is then for AAy = 1.2. Overall, the same diffi culty remains, how- 
ever. F igure[Halso shows the best-fit SED for Br uzual & Chariot! 
(2003) templates with the SMC law. Although the steepest SED 
between the optical and the H-band, it falls short in flux in the J 
band. For completeness we have also examined templates from 
the synthesis models of Marast on et alJ J2006). As expected, 
these models do not yield significantly different fits in the blue 
part of the rest-frame optical spectrum. 

As already mentioned by Lano rte et al.l d201 ll) we have also 
attempted to fit the SEDs with our m odels including nebular 
lines (see Sch aerer & de Barrosll20 09). Indeed, in this case the 
best-fit is found at z ~ 1.5 such that the [Om] AA 4959,5007 lines, 
and Hj6 boost somewhat the J-band flux, and Ha the H-band to 
a lesser extent, contributing thus to the flux decrement between 
J and Y. However, this solution requires also a very large atten- 
uation (Ay ~ 4.0 for the SMC law) to reproduce the steeply ris- 
ing SED towards longer wavelengths Q. Although to the best of 
our knowledge objects with such red SEDs and strong emission 
lines are not known, this extreme explanation should be easy to 
test with spectroscopic observations. 



4 Our model a ssumes identical attenuation for the continuum and 
nebular lines, as in ISchaerer & de Barrosl £2009) 



Finally, could composite populations not taken into account 
by our models help to explain the observed SED? Certainly 
the theoretical SED models may suffer from this simplification. 
However, we do not see how this could help to resolve the prob- 
lem with the large observed spectral break, since a superposition 
of individual simple stellar populations (not capable of repro- 
ducing this observation) can only average out spectral features. 
We conclude that we have no convincing explanation for the ob- 
served rapid drop of the observed SED of our two sources. 

5. Discussion 

5.1. Comparison with other objects in the literature 

How do our sources compare with other known galaxies and 
what is their nature? 

By design our sources are near-IR selected, optical drop- 
out sources, i.e. sources with a very red color between the J 
and z band and/or between Y and J . Our sources can there- 
fore be compared to those selected by Cap ak et alJ d201 ll) from 
the COSMOS survey. From their Fig. 13 we note that with 
(I-J) > 5.4 and 4.4 and (z-J) > 3.9 and 3.1 for zl and Y5, 
both sources show extreme (very red) optical to near-IR col- 
ors, when compared to other low redshift galaxies with red (J- 
z) colors, zl and Y5 are also similar to the z-dropout galaxy 
HUDF-J2 identified bv lMobasher et al l d2005l) as a z « 6.5 post- 
starburst galaxy candidate, but later shown to be most likely at 
z ~ 1 . 8-2.5 interloper by ISchaerer et ail (120071) : iDunlop et ail 
( 2007]) ; ICharv etakl d2007l) . Although similar in several respects, 
HUDF-J2 shows a more monotonously rising SED between z, J, 
and H than our objects exhibiting a "sharper" break. The col- 
ors and fluxes of our sources are also very similar to those of 
the lensed optical drop-out galaxy #2 found behind the clus- 
ter Abell 1835, identifi ed with the z = 2.93 sub-mm galaxy 
SMMJ 14009+0252 (see ISchaerer eTaIll2007t IWeifi et alJ 120091 
and references therein). 

Our sources are obviously also characterised by a large IR to 
optical flux ratio, a criterium which has been used by various au- 
thors. For example, Rodighiero et al. (2007) have studied IRAC 
3.6 fim selected sources undetected in deep optical HST images. 
The (K-3.6) and (z-3.6) colors of zl and Y5 are comparable to 
their sources; the main difference s eems to be in (H - K), w here 
our sources are bluer than those of [Rodighiero et al. (20Q3). zl 
and Y5 appear t o be related to the very dusty z ~ 2-3 so urces 
from this study. iFiore et al.l d2008l) and iDev et al.l (120081) have 
examined 24 /mi selected sources with very red colors between 
24 urn and the R band (S (24)/ S(R) Z 1000). They concluded 
that the bulk of these sources are very luminous strongly dust- 
obscured galaxies (referred to as DOGs) at z ~ 2, powered by 
AGN and/or by starbursts. Using the I band as a proxy for R, 
we obtain a flux ratio S(24)/S(R) ^ 20000 for zl, an extremely 
high flux ratio compared to the other samples. From the SED 
of Y5 (cf. Fig. [3), we also expect this galaxy to show a high 
24 /mi to optica l ratio. According to the source density from 
IDev et al.l (120081) we would have expected ~ 4 strongly dust- 
obscured galaxies down to ~ 0.3 mJy at 24 fim in our 45 arcmin 2 
field. The depth of our MIPS observations being similar, this 
value is compar able to our s ource density, although our selec- 
tion is different. iPope et aTl (120081) have also compared DOGs 
and sub-mm galaxies (SMGs) in the same R-K-24 /im colors, 
showing that ~ 30% of SMGs satisfy the DOG criteria, the re- 
mainder showing less extreme (i.e. bluer) colors. This confirms 
that the SEDs of our sources are comparable to a subset of SMGs 
with the most extreme optical to IR/sub-mm colors, as already 
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seen above (Sect.©. Our optical data is not deep enough to as- 
certain whether zl and Y5 fulfill the usual criteria for extremely 
red objects (EROs), (R - K) ^ 5.6 in Vega magnitudes, at least 
as estimated from (I-Ks). 

Among known sub-mm galaxies, one source, GN10 or 
GOODS 850-5, stands out as having particularly extreme 
IR/sub-mm to visible/near-IR properties, similar to our two 
galaxies. Indeed, this source is undetected down to ~ 0.01/Jy 
(1 cf) in the visible, shows fluxes of ~ 1 - 5/^Jy in the IRAC 
bands (3.6-8 nm) , and peaks at ~ 10 -20 mJy around 1 mm 
(IWang et al l 12004. 2009: iDaddi et al.l 12009-1) . quite comparable 
to zl and Y5. How ever, GN10 remains undetected even at JHK 
dWang et al.l2009l) . which can be explained by its hi gher redshift 
(z ~ 4 ), recently confirmed from CO spectroscopy (Dadd i~et"a"fl 
2009). The observed spectral break of GN 10 found between 
3.6 and 2.2 /urn and other considerations dWang et al.l [2009; 
Daddi et al. 2009) suggest a very high attenuation of Ay ^ 4.5-5 
for this source, or at least for the star-forming part of it, if host- 
in g multiple compone nts. If we assume constant star-formation 
as Daddi et al .1 (120091) for their SED modeling, we would infer 
Ay ~ 3 (7.8) mag for zl (Y5). Comparing the infrared-derived 
SFR with the upper limits in the rest-frame UV domain, we can 
also estimate the attenuation of our sources. Adopting the I-band 
flu x as a cons t raint f or the UV flux at ~ 2300-2500 A and using 
the lKennicuttl ([l998) calibration, we obtain Ay t A mag for both 
sources. 

Sub-mm galaxies are also known to exhibit very strong atten- 
uation. For example, source s with Balmer decrement measure- 
ments indicate Ay ~ 1-3 ( Swin bank et al. 20041: iTakata et al.l 
2006), and from SED fits (Swinb ank et al.ll2004l) es timate A v = 
3.0 + 1.0 for their sample. IWardlow et all (1201 ll) find Ay = 
2.6+0.2 from the median SED of sub-mm galaxies, but more ex- 
treme atte nuations are foun d within the sub-mm galaxy samples 
(see e.g. Smail et al. 2004]). Also, the e xtremely red obj ect, sub- 
mm detected galaxy HR10 studied by iDev et al.l d 19991) shows 
Ay ~ 4.5, as inferred from comparison of the IR and Ha star- 
formation rate. In short, although higher than the typical/median 
value of sub-mm galaxies, the attenuation of the sources zl and 
Y5 is similar to that of some sub-mm galaxies, such as GN10 at 
Z - 4.04 and others at lower z. Our sources are also somewhat 
fainter, both in the rest fra me near-IR and in the IR than the 
typical sub-mm galaxies (cf. IWardlow et al.ll20TTl) . Finally, our 
sources stand out by their large spectral break, which - to the 
best of our knowledge - is unusual among intermediate redshift 
sources. 



5.2. The other high-z candidates of the survey 

Based on our FIR detections we have identified two poten- 
tial intertor^e^samong the ten high-z candidates discovered 
by [Caporte et aT] d201 ll) . Most of the other candidates are in 
crowded regions where several sources emit in the FIR and are 
blended with each other, making any FIR measurement impos- 
sible. Two other candidates only seem to be clean from any con- 
tamination in the MIPS, PACS and SPIRE maps, namely Y3 and 
Y4. They remain undetected in all the bands. However, the FIR 
upper limits obtained do not allow us to discriminate between 
low and high redshift. On the other hand we can rule out that 
these sources are as extreme as zl and Y5 in their IR/sub-mm 
to near-IR flux ratio, since they should otherwise clearly be de- 
tected in our Herschel images. 



6. Conclusions 

Analyzing the FIR SED of two high redshift dropout candidates 
we find that both galaxies are likely at z ~ 2 rather than z > 7. 
From the FIR point of view alone, both galaxies could be similar 
to ULIRGs or SMGs which are common at z ~ 2. At z > 7 the 
SEDs would imply extreme dust temperatures and luminosities. 
Fitting the global SEDs considering all the data available from 
visible to submm we estimate z ~ 1.6-2.5. 

However, the optical/NIR part of both objects remains diffi- 
cult to understand if at z ~ 2. They show a very strong and well 
defined spectral break (presumably the Balmer break), unusual 
among intermediate redshift sources. The source zl is extreme 
for both IR/visible and near-IR/optical colors. The source Y5 has 
a somewhat smaller drop between the near-IR/optical domain. 

We have examined several possible explanations for the ex- 
treme colors of these galaxies but none of them is entirely satis- 
factory. More observations are required to understand their na- 
ture. Once a spectroscopic confirmation of their redshift is ob- 
tained we will be able to create new SED te mplates. 

Other extreme sources (e.g. GN10 IWang et all 120041: 



es ( e.g. 

iLutz et all I2005L IWang et al.1 120091: IDaddi et al.l 120091) were 
found from MIPS observations and submm searches. There 
could therefore be two complementary paths leading to similar, 
extreme sources. 

Although spectroscopic confirmation is still required and all 
possible interlopers may not be detected in the FIR, this work 
shows that FIR observations can be very helpful to constrain 
the contamination of high-z dropout searches by lower redshift 
galaxies. 
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